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Abstract. 

It has for long been an article of faith among astrophysicists 
that black hole spin energy is responsible for powering the 
relativistic jets seen in accreting black holes. Two recent 
advances have strengthened the case. First, numerical 
general relativistic magnetohydrodynamic simulations of 
accreting spinning black holes show that relativistic jets form 
spontaneously. In at least some cases, there is unambiguous 
evidence that much of the jet energy comes from the black 
hole, not the disk. Second, spin parameters of a number of 
accreting stellar-mass black holes have been measured. For 
ballistic jets from these systems, it is found that the radio 
luminosity of the jet correlates with the spin of the black hole. 
This suggests a causal relationship between black hole spin and 
jet power, presumably due to a generalized Penrose process. 



1. Introduction 

Relativistic jets are a common feature of accreting black holes (BHs) . They 
are found in both stellar-mass BHs and supermassive BHs, and are often 
very poviferful. Understanding how jets form and where they obtain their 
enormous power is an active area of research in astrophysics. 

In seminal work, Penrose [27] showed that a spinning BH has free 
energy that is, in principle, available to be tapped. This led to the 



popular idea that the energy source behind relativistic jets might be 
the rotational energy of the accreting BH. A number of astrophysical 
scenarios have been described in which magnetic fields enable this process 
[SSlElEllIIlinilTnillTllSllMlEl]. Field lines are kept confined around the 
BH by an accretion disk, and the rotation of space-time near the BH twists 
these lines into helical magnetic springs which expand under their own 
pressure, accelerating any attached plasma. Energy is thereby extracted 
from the spinning BH and is transported out along the magnetic field, 
making a relativistic jet. Although this mechanism requires accretion of 
magnetized fluid and is thus not the same as Penrose's original proposal, 
we will still refer to it as the "generalized Penrose process" since ultimately 
the energy comes from the spin of the BH. 

It is not easy to prove that the generalized Penrose process is necessarily 
in operation in a given jet. The reason is that jets are always associated 
with accretion disks, and the accretion process itself releases gravitational 
energy, some of which might flow into the jet. Let us define a jet efficiency 
factor 77jot, 

where (ijct) is the time-average (kinetic and electromagnetic) luminosity 
flowing out through the jet and (Af(rH))c^ is the time-average rate at 
which rest-mass energy flows in through the BH horizon. Many jets, both 
those observed and those seen in computer simulations, have values of ?7jct 
quite a bit less than unity. With such a modest efficiency, the jet power 
could easily come from the accretion disk [3J [HI [13] • 

The situation has improved considerably in the last couple of years. 
As we show in ^J21 numerical simulations have now been carried out where 
it can be demonstrated beyond reasonable doubt that the simulated jet 
obtains power directly from the BH spin energy. Furthermore, as we 
discuss in ^\ the first observational evidence for a possible correlation 
between jet power and BH spin has finally been obtained. The correlation 
appears to favor a Penrose-like process being the energy source of jets. 

^ Penrose considered a simple model in which particles on negative energy orbits 
fall into a spinning BH. jSTJ extended the analysis to discrete particle accretion in 
the presence of a magnetic field, which introduces additional interesting effects. 
We do not discuss these particle-based mechanisms, but focus purely on fluid 
dynamical processes within the magnetohydrodynamic (MHD) approximation. 
We also do not discuss an ongoing controversy on whether or not different 
mechanisms based on magnetized fluids differ from one another [12j . 
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2. Computer Simulations of Black Hole Accretion and Jets 

For the last decade or so, it has been possible to simulate numerically the 
dynamics of MHD accretion flows in the fixed Kerr metric of a spinning 
BH. The dynamics of the magnetized fluid are described using the general 
relativistic MHD (GRMHD) equations in a fixed space-time, and the 
simulations are carried out in 3D in order to capture the magnetorotational 
instability (MRI), the agency that drives accretion J^. Radiation is usually 
ignored, but this is not considered a problem since jets are usually found 
in systems with geometrically thick accretion disks, which are radiatively 
inefficient. We describe here one set of numerical experiments [Ml 135 
which have been run using the GRMHD code HARM [8] and which 
are particularly relevant for understanding the connection between the 
generalized Penrose process and jets. 

As is standard, the numerical simulations are initialized with an 
equilibrium gas torus orbiting in the equatorial plane of a spinning BH. 
The torus is initially embedded with a weak seed magnetic field, as shown 
in panel (a) of Figure [TJ Once the simulation begins, the magnetic field 
strength grows as a result of the MRI [T]. This leads to MHD turbulence, 
which in turn drives accretion of mass and magnetic field into the BH. The 
mass accretion rate is given by 

M{r) = -11 pu'-dAe,^, (2) 
J Je,ip 

where the integration is over all angles on a sphere of radius r (Boyer- 
Lindquist or Kerr-Schild coordinates), dAg^^ — ^J—gdOdip is the surface 
area element, p is the density, is the contravariant radial component of 
4- velocity, and g is the determinant of the metric. The sign in equation ^ 
is chosen such that M > corresponds to mass infiow. One is usually 
interested in the mass accretion rate at the horizon, M(r — ru). In 
computing M(rH) from simulations, one waits until the system has reached 
approximate steady state. One then computes M(rH) over a sequence 
of many snapshots in time and then averages to eliminate turbulent 
fluctuations. This gives the time-average mass accretion rate {Mir-n)). 

Panels (b)-(d) in Figure[T]show the time evolution of the accretion flow 
and jet in a simulation with BH spin a* = a/M — 0.99, where M is the 
BH mass [35]. The steady accretion of magnetized fiuid causes magnetic 
field to accumulate in the inner regions near the BH. After a while, the 
field becomes so strong that it compresses the inner part of the otherwise 
geometrically thick accretion flow into a thin sheet (panel b). The effect is 
to obstruct the accretion flow (panels c and d), leading to what is known 
as a magnetically-arrested disk [33] [T^ or a magnetically choked accretion 
flow |20| . The strong field extracts BH spin energy and forms a powerful 
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Figure 1. Formation of a magnetically-arrested disk and ejection of 
powerful jets in a GRMHD simulation of magnetized accretion on to a 
rapidly spinning BH with a* = 0.99 [3S]. The top and bottom rows in 
panels (a)-(d) show a time sequence of equatorial and meridional slices 
through the accretion flow. Solid lines show magnetic field lines in the 
image plane, and color shows logp (red high, blue low). The simulation 
starts with an equilibrium torus embedded with a weak magnetic field 
(panel a). The weakly magnetized orbiting gas is unstable to the MRI, 
which causes gas and field to accrete. As large-scale magnetic flux 
accumulates at the center, a coherent bundle of field lines forms at the 
center, which threads the BH and has the configuration of bipolar funnels 
along the (vertical) BH rotation axis. These funnels contain strong field 
and low mass density (lower panels b, c, d). Helical twisting of the field 
lines as a result of dragging of frames causes a powerful outflow of energy 
through the funnels in the form of twin jets. The outflow efficiency rj 
(panel e), calculated as in equation ([5|), becomes greater than unity once 
the flow achieves quasi-steady state at time t > 5000rg/c. This is the key 
result of the simulation. Having a time-average 77 > 1 means that there 
is a net energy flow out of the BH, i.e., spin energy is extracted from the 
BH by the magnetized accretion flow. This constitutes a demonstration of 
the generalized Penrose process in the astrophysically relevant context of 
a magnetized accretion flow. 
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outflow. To understand the energetics, consider the rate of flow of energy, 



E{t) ^ jj TldAe,^, (3) 
where the stress-energy tensor of the magnetized fluid is 

T^' =[p + u, +P, + ^) u^u. + + 5^ ~ ^, (4) 

Ug and Pg are the internal energy and pressure of the gas, is the fluid- 
frame magnetic field 4- vector (see [8] for the definition), and iP' = b'^bfi 
is the square of the fluid-frame magnetic field strength. The sign of 
equation ([3]) is chosen such that E{r) > corresponds to energy inflow. 
Note that includes the infiow of rest mass energy via the term pu^ut . 

Let us define the efficiency with which the accreting BH produces 
outfiowing energy as 

where we have made rj dimensionless by normalizing the right-hand side by 
the time-average mass energy accretion rate. To understand the meaning 
of equation (O, consider the simple example of gas falling in radially from 
infinity, with no radiative or other energy losses along the way. In this 
case, we have E{rii) = Af(rH)c^, i.e., the gas carries an energy equal 
to its rest mass energy into the BH. Hence 77 = 0, as appropriate for this 
example. For a more realistic accretion flow, some energy is lost by the gas 
via radiation, winds and jets, and one generally expects the energy flowing 
into the BH to be less than the rest mass energy: E{r}i) < M{rii)c^. This 
will result in an efficiency 77 > 0, where 77 measures the ratio of the energy 
returned to infinity, M(rH)c^ — -E(rH), to the energetic price paid in the 
form of rest mass energy fiowing into the BH, M{r}i)c^. 

Usually, E{rii) is positive, i.e., there is a net flow of energy into the 
BH through the horizon, and < 1. However, there is no theorem that 
requires this. Penrose's [2T great insight was to realize that it is possible 
to have E^rn) < (net outward energy flow as measured at the horizon), 
and thus 77 > 1. In the context of an accretion flow, E{r}i) < means that, 
even though rest mass flows steadily into the BH, there is a net energy flow 
out of the BH. As a result, the gravitational mass of the BH decreases with 
time. It is the energy associated with this decreasing mass that enables 77 
to exceed unity. Of course, as the BH loses gravitational mass, it also loses 
angular momentum and spins down. This can be verified by considering 
the angular momentum fiux at the horizon, J^rn), which may be computed 
as in equation but with T[ replaced by (e.g., |2S])- 
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Figure 2. Time-average outflow cfBciency 77 versus BH spin parameter 
a* for a sequence of GRMHD simulations of non-radiative BH accretion 
flows [34]. The efficiency exceeds unity for a* > 0.9. Negative values of a* 
correspond to the accretion flow counter-rotating with respect to the BH. 

Returning to the simulation under consideration, Figure [TJe) shows 
the outflow efficiency 77 as a function of time. It is seen that the 
average efficiency exceeds unity once the flow achieves steady state at 
time t > 5000rg/c, where Vg = GM/c^. The outflow thus carries away 
more energy than the entire rest mass energy brought in by the accretion 
flow. This is an unambiguous demonstration of the generalized Penrose 
process in the astrophysically plausible setting of a magnetized accretion 
flow on to a spinning BH. Of course, it is not obvious that the energy 
necessarily flows out in a collimated relativistic jet, since the quantity 77 is 
defined via global integrals (eqs. [51 [5]) which do not specify exactly where 
the outflowing energy ends up. A more detailed analysis reveals that the 
bulk of the energy does indeed go into a relativistic jet, while about 10% 
goes into a quasi-relativistic wind [34] . 

Figure [T] corresponds to an extreme example, viz., a very rapidly 
spinning BH with a* = 0.99. Figure [2] shows results from a parameter 
study that investigated the effect of varying a*. It is seen that the time- 
average rj increases steeply with increasing a*. For an accretion flow that 
corotates with the BH, the power going into the jet can be well-fit with a 
power-law dependence, 

?7jct«0.65a2(l + 0.85a2). (6) 

This approximation remains accurate to within 15% for 0.3 < a, < 1. 
For low spins, the net efficiency derived from the simulations is greater 
than that predicted by equation ((6|). For example, as Figure [2] shows, the 
simulation gives a non-zero value of rj for a* = 0, which is inconsistent 
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with equation ([5]). This is because, for a* — 0, all the outflow energy 
comes directly from the accretion how, most of which goes into a wind. 
Nothing comes from the BH, whereas equation ^ refers specifically to the 
efficiency T^jot associated with jet power from the BH. With increasing BH 
spin, both the disk and the hole contribute to energy outflow, with the 
latter becoming more and more dominant. For spin values a» > 0.9, the 
BH's contribution is so large that the net efficiency exceeds unity. 

Before leaving this topic, we note that other numerical simulations have 
used geometrically thicker accretion configurations than the one shown in 
Figure [T] and find even larger values of t] [33| [20] . 

3. Empirical Evidence for the Generalized Penrose Process 

As discussed in ^ there is definite evidence from computer simulations 
that the generalized Penrose process is feasible, and even quite plausible, 
with magnetized accretion flows. We discuss here recent progress on the 
observational front. In tj3.II we briefly summarize efforts to measure spin 
parameters of astrophysical BHs. Then in ^3.21 we discuss a correlation 
that has been found between jet power and BH spin. Finally in ^ 33. 31 
we explain why we think the observational evidence favors a Penrose-like 
process rather than disk power. 

3.1. Spin Parameters of Stellar-Mass Black Holes 

In 1989, the flrst practical approach to measuring black hole spin was 
suggested [3, viz., modeling the relativistically-broadened Fe K emission 
line emitted from the inner regions of an accretion disk. The first 
compelling observation of such a line was reported six years later [32] . 
Presently, the spins of more than a dozen black holes have been estimated 
by modeling the Fe K line (see [53] for a recent review). 

In 1997, a second approach to measuring black hole spin, the 
"continuum-fitting method," was proposed ^38 . In this method, one fits 
the thermal continuum spectrum of a black hole's accretion disk to the 
relativistic model of Novikov & Thorne [25]. One then identifies the 
inner edge of the modeled disk with the radius i?isco of the innermost 
stable circular orbit (ISCO) in the space-time metric. S ince -Risco varies 
monotonically with respect to the dimensionless BH spin parameter a* (see 
Fig. [3]) , a measurement of the former immediately provides an estimate of 
the latter. 

In 2006, the continuum-fitting method was employed to estimate the 
spins of three stellar-mass BHs [301 [IS] • Seven additional spins have since 
been measured. Table [T] lists the masses and spins of these ten BHs. 
Readers are referred to a recent review ■ 14; for details of the continuum- 
fitting method and uncertainties in the derived spin estimates. 
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Figure 3. Radius of the ISCO i?isco and of the horizon i?H in units of 
GMjc} plotted as a function of the black hole spin parameter a, . Negative 
values of a* correspond to retrograde orbits. Note that i?isco decreases 
monotonically from ^GMj (? for a retrograde orbit around a maximally 
spinning black hole, to QGAI/c^ for a non-spinning black hole, to GM/(? 
for a prograde orbit around a maximally spinning black hole. 

The continuum- fitting method is simple and demonstrably robust. It 
does not make many assumptions; those few it makes have nearly all 
been tested and shown to be valid (see [16j[T4] for details). A significant 
limitation of the method is that it is only readily applicable to stellar- 
mass BHs. For such BHs, however, we would argue that it is the method 
of choice. The Fe K method can be applied to both stellar-mass and 
supermassive BHs. For the latter, it is the only method currently available. 

3.2. Correlation Between Black Hole Spin and Jet Radio Power 

The 10 stellar-mass BHs in Table [1] are divided into two classes: 
"persistent" sources, which are perennially bright in X-rays at a relatively 
constant level, and "transient" sources, which have extremely large 
amplitude outbursts. During outburst, the transient sources generally 
reach close to the Eddington luminosity limit (see [31] for a quantitative 
discussion of this point). Close to the peak, these systems eject blobs 
of plasma that move ballistically outward at relativistic speeds (Lorentz 
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Table 1. The spins and masses of ten stellar- mass black holes [H] . 



System 


a. 


M/Mq 


Persistent 


Cygnus X-1 


> 0.95 


14.8 ± 1.0 


LMC X-1 




10.9 ± 1.4 


M33 X-7 


0.84 ±0.05 


15.65 ± 1.45 


Transient 


GRS 1915-H05 


> 0.95 


10.1 ±0.6 


4U 1543-47 


0.80 ±0.10 


9.4 ± 1.0 


GRO J1655-40 


0.70 ±0.10 


6.3 ±0.5 


XTE J 1550-564 


0.34 ± 0.24 


9.1 ±0.6 


H1743-322 


0.2 ±0.3 


~ 8 


LMC X-3 


< 0.3 


7.6 ±1.6 


A0620-00 


0.12 ±0.19 


6.6 ±0.25 



factor r > 2). These ballistic jets are often visible in radio and sometimes 
in X-rays out to distances of order a parsec from the BH, i.e., to distances 
> IQ^^GM / (? . Because ballistic jets resemble the kiloparsec-scale jets seen 
in quasars, stellar-mass BHs that produce them are called microquasars 
[22]. 

On general principles, one expects the jet luminosity Ljct (the power 
flowing out in the kinetic and electromagnetic energy) to depend on the BH 
mass M, its spin a*, and the mass accretion rate M (plus perhaps other 
qualitative factors such as the topology of the magnetic field [21 HH])- If one 
wishes to investigate the dependence of jet luminosity on a*, one needs first 
to eliminate the other two variables. Ballistic jets from transient stellar- 
mass BHs are very well-suited for this purpose. First, the BH masses are 
similar to better than a factor of two (see Table [l}. Second, all these 
sources have similar accretion rates, close to the Eddington limit, at the 
time they eject their ballistic jets |31j . This leaves a, as the only remaining 
variable. 

Narayan & McClintock [53] considered the peak radio luminosities 
of baUistic jet blobs in four transient BHs, A0620-00, XTE J1550-564, 
GRO J1655-40, GRS 1915±105, and showed that they correlate well with 
the corresponding black hole spins measured via the continuum-fitting 
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Figure 4. (a) Plot of jet power estimated from 5 GHz radio flux at light 
curve maximum, versus black hole spin measured via the continuum-fitting 
method, for five transient stellar-mass BHs [M] |3T] . The dashed line has 
slope fixed to 2 (see eq. ^ and is not a fit. (b) Plot of jet power versus 
-Risco / [GM/ c^). Here jet power has been corrected for beaming assuming 
jet Lorentz factor F = 2 (filled circles) or F = 5 (open circles). The two 
solid lines correspond to fits of a relation of the form, "Jet Power" oc $7^ , 
where fin is the angular frequency of the black hole horizon [3T]. Note 
that jet power varies by a factor > 1000 among the five objects shown. 

methocd. Later, "HT included a fifth BH, H1743-322, whose spin had been 
just measured. Figure shows the inferred ballistic jet luminosities of 
these five objects plotted versus black hole spin. The quantity "Jet Power" 
along the vertical axis refers to the scaled maximum radio luminosity 
{vSu)D'^ /M , where v (= 5 GHz) is the radio frequency at which the 
measurements are made, is the radio flux density in janskys at this 
frequency at the peak of the ballistic jet radio light curve, D is the 
distance in kiloparsecs, and M is the black hole mass in solar units. There 
is unmistakeable evidence for a correlation between the observationally 
determined quantity Jet Power and a*. Although there are only five data 
points, note that Jet Power varies by nearly three orders of magnitude as 
the spin parameter varies from « 0.1 — 1. 

The very unequal horizontal errorbars in Figure |4ji are a feature of 
the continuum- fitting method of measuring a, . Recall that the method in 
effect measures -Risco a-nd then deduces the value of a* using the mapping 

^ In the case of a fifth transient BH, 4U1543-47, radio observations did not 
include the peak of the light curve, so one could only deduce a lower limit to the 
jet power, which is shown as an open circle in Figure 2^. 
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shown in Figure [31 Since the mapping is non-hnear, especiahy as a^, — 1, 
comparable errors in i?isco correspond to very different uncertainties 
in a*. In addition, the use of log a* along the horizontal axis tends to 
stretch errorbars excessively for low spin values. Figure IHd, based on [31], 
illustrates these effects. Here the horizontal axis tracks logi?isco rather 
than log a*, and the horizontal errorbars are therefore more nearly equal. 
The key point is, regardless of how one plots the data, the correlation 
between Jet Power and black hole spin appears to be strong. 

3.3. Why Generalized Penrose Process? 

Assuming the correlation shown in Figure[l|is real, there are two immediate 
implications: (i) Ballistic jets in stellar-mass BHs are highly sensitive to 
the spins of their underlying BHs. (ii) Spin estimates of stellar-mass BHs 
obtained via the continuum-fitting method are sufficiently reliable to reveal 
this long-sought connection between relativistic jets and BH spin. 

With respect to (i), the mere existence of a correlation does not 
necessarily imply that the generalized Penrose process is at work. We 
know that the accretion disk itself is capable of producing a jet-like outflow 
[H [HI [13]. Furthermore, the gravitational potential well into which an 
accretion disk falls becomes deeper with increasing BH spin, since the 
inner radius of the disk i?isco becomes smaller (Fig. [3]). Therefore, a 
disk-driven jet is likely to become more powerful with increasing spin. 
Could this be the reason for the correlation between Jet Power and spin 
seen in Figure [D? We consider it unlikely. The radiative efficiency jydisk of a 
Novikov-Thorne thin accretion disk increases only modestly with spin; for 
the spins of the five objects shown in Figure [H rydisk = 0.061, 0.069, 0.072, 
0.10 and 0.19, respectively, varying by only a factor of three. Of course, 
there is no reason why the power of a disk-driven jet should necessarily 
scale like r^disk- Nevertheless, the fact that Tydisk shows only a factor of 
three variation makes it implausible that a disk-powered jet could have its 
radio luminosity vary by three orders of magnitude. 

In contrast, any mechanism that taps directly into the BH spin energy 
via some kind of generalized Penrose process can easily account for 
the observed variation in Jet Power. Analytical models of magnetized 
accretion predict that the jet efficiency factor should vary as Tyjot oc 
[ini [S] or 77jet oc [3S] , where fin is the angular frequency of the BH 
horizord. The dashed line in Figure [1^ corresponds to the former scaling, 
and the sohd lines in Figure [Hd to the latter scaling; equation which 
is obtained by fitting simulation results, is intermediate between the two. 
The observational data agree remarkably well with the predicted scalings, 
strongly suggesting that the generalized Penrose process is in operation. 

^ The two scalings agree for small values of a,, but differ as a, — s- 1. 
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We cannot tell whether the energy extraction in the observed systems 
is mediated specifically by magnetic fields (as in the simulations), since 
there is no way to observe what is going on near the BH at the place 
where the jet is initially launched. Where the ballistic jet blobs are finally 
observed they are clearly magnetized — it is what enables the charged 
particles to produce radiation via the synchrotron mechanism — but this 
is at distances IQ^^GM/c^. 

4. Summary 

In summary, the case for a generalized version of the Penrose process 
being the power source behind astrophysical jets has become significantly 
stronger in the last few years. Computer simulations have been very helpful 
in this regard since they enable one to study semi-realistic configurations 
of magnetized accretion flows and to explore quantitatively how mass, 
energy and angular momentum flow through the system. Recent computer 
experiments find unambiguous indications for energy extraction from 
spinning BHs via magnetic fields. Whether these simulated models 
describe real BHs in nature is not yet certain. However, completely 
independent observational data suggest a link between the spins of 
transient stellar-mass BHs and the energy output in ballistic jets ejected 
from these systems. The observationally measured quantity "Jet Power" in 
Figure m increases steeply with the measured BH spin, and the dependence 
is quite similar to that found both in simple analytical models |291 [5] 
and in simulations (Fig. [2]). Taking all the evidence into account, the 
authors believe that Penrose's seminal ideas on energy extraction from 
spinning BHs are relevant for the production of at least some categories of 
relativistic astrophysical jets. 

The authors thank J. F. Steiner for help with Fig. HJd. RN's work 
was supported in part by NASA grant NNX11AE16G. AT was supported 
by a Princeton Center for Theoretical Science fellowship and an XSEDE 
allocation TG-AST100040 on NICS Kraken and Nautilus and TACC 
Ranch. 
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